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Abstract 

The positively charged pions produced in proton-proton collisions at a beam mo- 
mentum of 1640 MeV/c were measured in the forward direction with a high resolu- 
tion magnetic spectrograph. The missing mass distribution shows the bound state 
(deuteron) clearly separated from the pn continuum. Despite the very good reso- 
lution, there is no evidence for any significant production of the pn system in the 
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spin-singlet state. However, the a{pp -K'^pn) / a{pp n^d) cross section ratio is 
about twice as large as that predicted from S-wave final-state-interaction theory 
and it is suggested that this is due to D-state effects in the pn system. 
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There is a very extensive literature on the pp Ti^d reaction and many 
detailed analyses have been made [1], but much less is known about the pro- 
duction of the continuum in the pp — > vr+pn case. Data covering low exci- 
tation energies generally show the strong 5"- wave final-state-interaction (fsi) 
peak corresponding to the pn spin-triplet which has, as a characteristic energy 
scale, the binding energy of the deuteron {Bt = 2.22 MeV). However, the en- 
ergy resolution is generally insufficient to identify the analogous spin-singlet 
fsi peak, for which the corresponding energy scale is only Bg = 0.07 MeV [2]. 
Indirect evidence suggests that spin-singlet production is much weaker than 
that of spin-triplet for medium energy proton beams [3] , and this is confirmed 
by data from the isospin-related pp — > n^pp reaction, though these arc lim- 
ited in incident momentum or energy resolution [4]. Such weak spin-singlet 
production accords well with theory, because the influence of the A-isobar is 
minimal there. 

A useful way of trying to extract the spin-singlet contribution is through the 
comparison of the overall strengths of the cross sections for pn and deuteron 
final states. Using final-state-interaction theory, Faldt and Wilkin derived the 
extrapolation theorem which relates the normalisations of the wave functions 
for 5'-wave bound and scattering states [5] . This has been exploited to predict 
the double-differential centre-of-mass (cm) cross section for the S-wave spin- 
triplet component in pp — > ir'^pn in terms of the cross section for pp n^d [6]: 



dQdx'"^"^ p(-l) 27r(a; + 1) dQ 



Here x denotes the excitation energy e in the np system in units oi Bt, x ~ 
e/Bt, and p{x) and p(— 1) are the pion cm momenta for the pn continuum or 
deuteron respectively. 

In the derivation of Eq. (1) it is assumed [6] that the pion production operator 
is of short range and that x is not too large, so that the pn P- waves contribute 
little. Most critical though is the neglect of channel coupling through the pn 
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tensor force, so that the equation could only be valid provided that the D-state 
effects are small in the production of both the bound state and continuum. 

The fsi peak arises from the ^/x/ (x + l) factor in Eq. (1) and there should be 
an analogous spin-singlet enhancement, where the deuteron binding energy Bt 
is replaced by the energy Bs of the virtual state in the S — 0, T — 1 system. 
At low excitation energies one therefore expects that 

^{PP - vr WJ = e ^{PP - n {pn},) , (2) 

where we use the factor ^ to quantify the ratio of spin-singlet to spin-triplet 
production. 

Since the best resolution in excitation energy so far achieved was typically 
a — 350 keV [7], any singlet peak would have been smeared significantly in all 
published data. However, by estimating the 5- wave triplet contribution to the 
pp TT+pn cross section from Eq. (1) and subtracting it from the observed 
data, some measure for the singlet production could be obtained. In most ex- 
periments where only the tt"*" was detected, the limited resolution did not guard 
against some leakage of the deuteron peak into the continuum region [8,9,10]. 
On the other hand, detecting the tt"'" and proton in coincidence [11], while 
identifying well the continuum channel, loses the relative normalisation with 
the TT^d final state, which is so important in the implementation of Eq. (1). 
Therefore, in addition to the pion spectrum, Betsch et al. [12] measured co- 
incidences between pion and proton, but then had to rely on Monte Carlo 
simulations. For 600 MeV and below, the data seemed to confirm that the 
singlet contributed at most 10% of the cross section, though at 1 GeV a higher 
figure was likely [9]. 

Most of the uncertainties mentioned above could be minimised by measur- 
ing just the pion spectrum, but with high resolution. One could then identify 
clearly any singlet peak and also separate unambiguously the pp n^pn from 
the pp — > TT^d reaction. This was our primary goal when planning a new exper- 
iment. Pions were observed near zero degrees with the 3Q2D spectrograph Big 
Karl [13] at the COSY accelerator in Jiilich. Their position and track direction 
in the focal plane were measured with two packs of multiwire drift chambers, 
each having six layers. The chambers were followed by scintillator hodoscopes 
that determined the time of flight over a distance of 3.5m. In order to optimise 
the momentum resolution, a liquid hydrogen target of only 2mm thickness was 
used with windows made of 1 /xm Mylar [14]. The beam was electron cooled at 
injection energy and, after acceleration, stochastically extracted. This resulted 
in an energy resolution of a = 97 keV for the deuteron peak. This was much 
better than that found in a test run without beam cooling and, in particular, 
the background was considerably reduced. 
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Fig. 1. The results from the present experiment (histogram) compared with the 
prediction (curve) of the -S-wave fsi theory of Eq. (1) [6]. 

The results of our experiment are shown in Fig. 1 as function of the excitation 
energy in the pn system. Though corrections for acceptance etc. have been 
included, these in fact vary slowly with e for energies below 20MeV. Noting the 
logarithmic scale in the figure, it is clear that there is an excellent distinction 
between the pp — > Tr'^pn from the pp — > 7r+d reactions. Since the luminosity 
and detection efficiencies largely cancel out between them, this means that 
we have a very good determination of the relative cross sections for ir'^d and 
TT+pn final states. 



Also shown in Fig. 1 is the prediction of the continuum production from the 
(S-wave fsi theory of Eq. (1), where we have assumed a constant background of 
30 counts per bin. Though the shape is largely right, it is too low in magnitude 
by a factor of 2.2 ±0.1 over the whole of the spectrum. It is interesting to note 
that, if our data are artificially degraded such that the resolution is the same 
as that achieved in the Leningrad experiment at the neighbouring energy of 
1 GeV ((7 fa 3 MeV) [9], the two sets of results overlap very well. However, 
the poor resolution allowed the authors of ref. [6] to ascribe the factor-of-two 
discrepancy to the production of spin-singlet final states. We can, however, 
check this hypothesis independently by studying the shape of the missing- 
mass spectrum. 

As is evident from Eq. (2), the cross section for producing a pn singlet state 
must show a sharp spike just above threshold and, due to our good resolu- 
tion, this prominent feature should remain even after convolution with this 
resolution. In Fig. 2 are shown the predictions of Eqs. (1,2) modified by the 
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Fig. 2. Comparison of the measured pn excitation energy spectrum on a linear scale 
with the prediction of Eqs. (1,2) for the shape of the singlet cross section. The error 
bars contain a tiny contribution from the uncertainty in the acceptance correction. 

inclusion by an extra factor of {l + s/Eg) to try to take into account deviations 
from the extrapolation theorem [5] . The value of Eg — 24 MeV is derived from 
the scattering length and effective range [15] though, by the point that this 
becomes significant, the S'-wave ansatz is dubious. There is no hint of any 
sharp needle in the data of Fig. 2 and, in fact, the shape of the cross section is 
completely compatible with pure spin-triplet production. Fits of Eq. (2) with 
free amounts of singlet and triplet show that ^ < 10~^ at the one standard 
deviation level, and this corresponds to a practically vanishing fraction if the 
singlet part. As a consequence, we must seek elsewhere for the factor-of-two 
discrepancy between our data and the results of Eq. (1). 

As has been stressed previously, the extrapolation theorem linking the bound 
and scattering wave functions is only valid if one can neglect completely D- 
state effects [5]. Though the D-state wave functions are suppressed at short 
distances by the centrifugal barrier, the ^-wave is also reduced in this region 
by the repulsive core. Thus the D-state might be significant for pion produc- 
tion despite the relatively small probability in the deuteron, especially if S-D 
interference terms are important. 

We consider a microscopic calculation of the actual three-body n^pn final 
state reaction to be beyond the scope of the present work. Nevertheless, to 
investigate the effects of the D-waye, at least semi-quantitatively, we have 
made estimations of the pp vr+rf differential cross section following the 
formalism described in ref. [16]. Using a standard deuteron wave function [17] 
with a normal D state, this reproduces well the experimental data [1]. The 
calculations have, however, been repeated with a reversed sign for the D-state 
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Fig. 3. D-state effects in tfie predicted excitation function [16] for the zero degree 
pp — >■ TT+d differential cross section. Tlie solid curve shows the results with the 
standard value [17], the broken curve with the reversed sign, and the dots with no 
D-state at all. 

amplitude and also with no £)-state at all. Now for kinematic reasons the pn 

D-state scattering wave function must vanish like £^ as £ — so that its sign 
should change when going from the bound state (deuteron) to the continuum 
pn pair [5]. One can therefore get an idea of the effect of the D-state in the 
continuum by using a deuteron wave function with the opposite sign for the 
D wave. 

The predictions for the forward cross section are shown in Fig. 3 as a function 
of the dimensionless pion cm momentum 7] = p/mTr+, the present experiment 
corresponding to = 2.6. The zero D-state calculation is approximately the 
average of the other two, showing that the effects are mainly due to S-D 
interference. At low energies the inclusion of the D state increases the cross 
section, while the converse is true at high energies. The exact position of the 
cross-over point, here predicted to be at ~ 1.6 (Tp 600 MeV), could 
be model dependent but we would certainly expect there to be a different 
influence of the £)-state on either side of the A peak. Given the theoretical 
uncertainties, the fact that the factor of 2.2 difference between the calculations 
with the changed sign of the D-state at rj = 2.6 coincides exactly with the 
discrepancy between the data and the ^-wave theory shown in Fig. 1 may 
be fortuitous. Close to or just below the resonance one would expect smaller 
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deviations from the extrapolation theorem associated with the D-state, and 
this certainly seems to be the case experimentally [6,7,8,10]. To quantify the 
deviations would require further high resolution runs which could identify 
clearly the singlet production from the shape of the spectrum. 

In summary, we have measured the missing mass spectrum from the pp —> 
7r~^X reaction in the forward direction. Despite the rather high beam momen- 
tum of 1640 MeV/c, the excellent resolution allowed the complete separation 
of the deuteron from pn continuum and also showed that the production of 
spin-singlet states was negligible at this momentum. Deviations from the re- 
sults of S'-wave fsi theory could be ascribed semi- quantitatively to the effects 
of the tensor force in the pn system and an extension of this to encompass the 
coupled S-D system would be of great help. It is also to be hoped that a full 
microscopic calculation of the three-body Tr'^pn final state production will be 
undertaken to complement the two-body results quoted here [16]. This might 
then confirm our hypothesis of the great influence of the deuteron I?-state in 
pion production above the A resonance. 
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